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INTRODUCTION 


Recent  advances  In  solid  propellant  rocket  motor  combustion  stability 
analysis  now  make  It  possible  to  calculate  the  stability  of  a proposed 
motor  design.1'2  Each  of  the  three  solid  propellant  rocket  motor  designs 
to  be  used  on  Space  Shuttle  flights  for  the  National  Aeronautics  and  Space 
Administration  (NASA)  is  to  have  stability  prediction  techniques  applied 
to  it  as  early  in  the  design  phase  as  feasible.  A general  flow  chart 
outlining  the  principal  items  involved  in  a stability  analysis  is  shown  in 
Figure  1.  As  indicated  in  that  figure,  detailed  knowledge  of  the  internal 
acoustic  characteristics  of  a motor  is  a necessary  prerequisite  to  deter- 
mination of  the  motor's  stability.  This  report  is  concerned  with  the 
acoustic  characteristics  of  the  largest  of  the  three  Space  Shuttle  solid 
propellant  motors:  the  solid  propellant  booster  motor  (SRM) . 

A matter  of  concern  unique  to  the  Space  Shuttle  relates  to  the  fact 
that  the  shuttle  vehicle  involves  two  solid  propellant  booster  motors  in 
conjunction  with  three  liquid  fuel  engines,  all  of  which  provide  thrust 
during  the  launch  phase.  There  is  a possibility  that  a standing  acoustic 
wave  in  the  solid  propellant  booster  motor  could  mechanically  couple 
through  the  vehicle  structure  with  the  liquid  fuel  or  oxidizer  feed 
system  to  create  an  oscillation  in  the  rate  of  injection  into  the  liquid 
engines.  Previous  experience  with  liquid  fueled  rocket  engines  shows  that 
a situation  can  occur  in  which  the  liquid  injection  fluctuations  cause 
engine  thrust  perturbations  which  in  turn  reinforce  and  amplify  the  liquid 
injection  variations.  This  behavior  is  known  as  the  POGO  effect  and  it  is 
an  undesirable  feature  of  engine  operation. 

Thus,  knowledge  of  the  internal  acoustic  characteristics  of  the  SRM 
is  not  only  necessary  for  the  solid  propellant  motor  stability  analysis 
but  it  is  also  required  for  determining  the  effect  of  internal  standing 
acoustic  waves  on  structural  response  such  as  dynamic  loading  of  the 
propellant  grain,  flexures  of  the  motor  case,  and  force  perturbations  on 
the  nozzle  assembly.  In  addition,  knowledge  of  the  frequencies  likely  to 
be  generated  by  the  SRM  can  provide  useful  data  to  the  liquid  engine 
system  designers  who  can  apply  that  information  to  existing  techniques  for 
reducing  the  probability  of  having  a POGO  effect  in  the  shuttle. 


Chemical  Propulsion  Information  Agency.  "Acoustic  Stability  Charac- 
terization of  the  Trident  (C-4)  Motors,"  by  M.  W.  Beckstead,  et.al.,  11th 
JANNAF  Combustion  Meeting.  Silver  Spring,  Md.,  CPIA,  December  1974,  p. 
535.  (CPIA  Pub.  261,  Vol.  I,  publication  UNCLASSIFIED.) 

2 . "Computer  Programs  for  Solid  Rocket  Motor  Stability  Predic- 

tions," by  R.  L.  Lovine  and  R.  C.  Waugh.  12th  JASSAF  Combustion  Meeting. 
Silver  Spring,  Md.,  CPIA,  December  1975,  p.  1.  (CPIA  Pub.  273,  Vol.  II, 
publication  UNCLASSIFIED.) 
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Analysis  of  the  internal  acoustics  of  the  SRM  was  conducted  in  two 
phases.  The  initial  phase  involved  use  of  equations  for  acoustic  fre- 
quency derived  from  the  classical  wave  equation  along  with  the  assump- 
tion that  the  interior  motor  geometry  could  be  approximated  by  a right 
circular  cylinder  (details  of  these  calculations  and  results  are  con- 
tained in  Appendix  A).  However,  the  internal  geometry  of  the  SRM  departs 
significantly  in  some  respects  from  the  idealized  cylindrical  form 
assumed  in  the  approximate  (initial)  calculation  as  shown  in  Figure  2. 
The  second  phase  of  the  analysis  of  SRM  internal  acoustics,  which  is  the 
main  subject  of  this  report,  involved  the  use  of  a relatively  sophisti- 
cated computer  method  to  provide  more  accurate  predictions  of  frequency 
and  of  acoustic  wave  structure  than  could  be  obtained  by  the  use  of 
classical  acoustics.  In  addition  to  providing  accurate  acoustic  wave 
characteristics  for  a non-cylindrical  interior,  the  computer  method  of 
acoustic  analysis  is  an  integral  part  of  the  motor  stability  calcula- 
tion. 


INPUT  DATA  AND  METHOD  OF  CALCULATION 


Analysis  of  the  acoustic  characteristics  of  a complicated  motor 
geometry  such  as  the  SRM,  using  presently  available  techniques,  involves 
use  of  a finite  element  formulation  of  the  problem  which  is  solved  with 
the  aid  of  a large,  high-speed  digital  computer.  The  NASTRAN  program, 
originally  developed  for  NASA  to  solve  problems  in  structural  dynamics, 
provides  a well-established  finite  element  technique  which  has  been 

* adapted  to  solving  the  problem  of  determining  the  natural  standing 
acoustic  waves  in  cavities  which  deviate  from  the  geometry  of  an  ideal 
cylindrical  shape.3*4  Two  methods  using  the  NASTRAN  program  are  avail- 
able at  the  Naval  Weapons  Center  (NWC) . One  method  involves  a quasi- 
three-dimensional  (3D)  program  which  requires  that  the  central  region  of 
the  cavity  be  circular  in  cross-section,  that  the  central  region  com- 
prise most  of  the  cavity  volume,  that  the  symmetry  be  cyclic,  and  that 

slots  radiating  from  the  central  cavity  be  narrow  in  relation  to  the  . 

cavity  diameter.  The  other  method  developed  from  the  3D  method  by  the 
second  author,  solves  the  acoustics  problem  in  two-dimensions.  The  two-  i • 

dimensional  (2D)  method  does  not  have  the  narrow  slot  restriction  which 
is  contained  in  the  3D  method.  Both  methods  require  similar  input 
Information  which  includes:  cavity  geometry,  boundary  conditions,  and 
parameters  relating  to  properties  of  the  gas  filling  the  cavity.  For 

L‘  • . f 

National  Aeronautics  and  Space  Administration.  NASTRAN  User's 

• Manual  (Level  15).  NASA,  June  1972.  (Publication  UNCLASSIFIED.) 

4National  Aeronautics  and  Space  Administration.  NASTRAN  Theoretical 
Manual  (Level  IS).  NASA,  April  1972.  (Publication  UNCLASSIFIED.) 
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rocket  motors,  the  combustion  gas  properties  are  either  obtained  from 
data  provided  by  the  propellant  manufacturer  or  through  the  use  of  a 
computer  Propellant  Evaluation  Program  (PEP)  which  is  available  at  NWC. 
If  the  PEP  method  is  used,  as  in  the  case  of  the  SRM,  the  propellant 
formulation  and  the  motor  operating  pressure  are  required  as  inputs  to 
the  program. 

The  SRM  meets  the  internal  geometry  requirements  of  the  3D  NASTRAN 
method  for  determining  acoustic  characteristics.  Therefore,  since  both 
axial  and  transverse  acoustic  modes  were  of  interest,  it  was  determined 
that  the  proper  approach  was  to  use  the  3D  program. 

Solution  of  the  acoustic  characteristics  of  the  booster  motor 
requires  setting  up  a finite-element  grid  for  each  time  during  burn. 
Grids  for  three  internal  configurations  were  established  using  large 
scale  drawings  of  the  motor  at  zero  web  burn  (provided  by  Thiokol, 
Wasatch  Division)  and  drawings  of  the  propellant  surface  regression 
(furnished  by  Rockwell  International).  The  three  configurations'  repre- 
sented web  burns  of  0,  48,  and  86  cm  (0,  19,  and  34  inches),  respective- 
ly. The  grid  system  used  for  the  0-cm  web  burn,  shown  in  Figure  3,  is 
typical  of  those  used  for  the  other  internal  configurations.  Each  of 
the  grids  used  is  shown  in  detail  in  Appendix  B and  grid  coordinate  data 
for  each  of  the  three  configurations  are  provided  in  Appendix  C. 

An  Important  assumption  in  the  acoustic  analysis  is  that  the 
boundaries  of  the  gas-filled  interior  of  the  motor  are  treated  as  rigid 
walls.  Thus,  the  acoustic  mode  program  did  not  allow  for  transfer  of 
energy  from  the  gas  oscillations  to  the  propellant  grain  or  to  the  motor 
case.  Another  assumption  used  in  the  present  analysis  is  that  the  speed 
of  sound  in  the  gas  is  uniform  throughout  the  cavity.  In  addition,  no 
allowance  is  made  in  the  program  for  mean  gas  flow.  However,  the  pro- 
gram does  allow  for  the  assumption  to  be  made  of  a closed  or  an  open 
nozzle  throat.  Both  assumptions  have  been  used  in  determining  the  SRM 
acoustic  characteristics  as  explained  in  the  following  paragraphs. 

It  has  been  the  custom,  in  assessing  rocket  motor  acoustic  char- 
acteristics, to  assume  a closed  throat  condition.  The  assumption  is  a 
convenience  in  chat  it  simplifies  the  analysis,  particularly  when  per- 
forming rapid  hand  calculations  (see  Appendix  A).  Furthermore,  the 
closed  throat  assumption  is  a reasonable  one  for  the  many  rocket  motors 
which  have  a small  ratio  of  nozzle  throat  area  to  propellant  gas  port 
area.  This  area  ratio,  usually  termed  "J",  is  not  small,  however,  in 
the  SRM:  it  is  initially  0.716  at  Ignition  and  drops  to  a value  of 
0.381  at  burnout. 

The  issue  of  how  to  treat  the  nozzle  throat  acoustically  for  the 
SRM  is  not  clear  at  the  present  time.  One  authority  suggests  that 
converging  gas  flow  in  the  nozzle  entry  region  is  a possible  source  of 
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ncoustic  reflection  and  that  It  is  therefore  plausible  to  treat  the 
nozzle  throat  as  an  acoustic  reflector  (or  as  a rigid  surface).6  On 
the  other  hand,  experimental  cold  gas  flow  data  in  small  scale  model 
rocket  motors  indicates  that  the  axial  acoustic  wave  characteristics  of 
the  motor  are  affected  by  J in  a way  that  suggests  the  nozzle  should  be 
treated  as  an  open  or  non-reflecting  area.6  A brief  discussion  of  the 
experimental  cold  flow  data  is  found  in  Appendix  D. 

An  additional  factor  which  is  known  to  have  an  etfect  on  the  in- 
ternal acoustics  of  the  SRM  is  the  flow  of  gas  in  the  motor  and  nozzle 
entry  area.  However,  information  regarding  gas  flow  in  the  motor  is  not 
included  in  the  present  acoustic  calculations  so  that  a quantitative 
assessment  of  its  effect  on  motor  acoustics  is  not  possible  at  this 
time.  In  regard  to  gas  flow,  its  importance  increases  as  J increases  so 
that  a motor  with  large  J will  have  two  related  phenomena  affecting  its 
internal  acoustics:  (1)  an  essentially  geometric  effect  caused  by 
presence  of  the  nozzle  throat  as  an  area  in  which  acoustic  waves  are 
"absorbed",  and  (2)  the  effect  of  a net  gas  flow  in  the  motor-  which 
introduces  an  "asymmetry"  into  the  propagation  of  acoustic  waves:  waves 
traveling  in  the  direction  of  flow  (toward  the  nozzle)  have  a greater 
velocity  than  waves  moving  upstream  (toward  the  head  end  of  the  motor). 

In  light  of  the  present  situation  regarding  treatment  of  the  acous- 
tic behavior  of  the  nozzle  it  seems  prudent  to  present  acoustic  data  for 
both  throat  conditions.  However,  as  noted  above  neither  set  of  data 
includes  the  effects  of  mean  gas  flow  on  acoustics. 


RESULTS 


The  3D  NASTRAN  acoustics  program  provides  a combination  of  tabu- 
lated and  graphic  output  data.  The  tabulated  data  consists  primarily  of 
an  acoustic  pressure  distribution  which  is  in  the  form  of  a normalized 
pressure  for  each  point  in  the  finite  element  grid  system.  For  each 
acoustic  wave  solution,  the  tabulated  pressures  are  expressed  as  frac- 
tions of  the  maximum  pressure  in  the  cavity  (which  is  assigned  a value 
of  unity).  The  graphic  output  provides  a plan  view  of  the  finite  element 


^Chemical  Propulsion  Information  Agency.  "Combustion  Instability 
in  Large  Solid  Rocket  Motors,"  by  F.  E.  C.  Culick  and  R.  N.  Kumar,  10th 
JANNAF  Combustion  Meeting.  Silver  Spring,  Md.,  CPIA,  December  1973,  p. 
45.  (CPIA  Pub.  243,  Vol.  I,  publication  UNCLASSIFIED.) 

®F.  G.  Buffum,  Jr.,  G.  L.  Dehority,  R.  0.  Slates,  and  E.  V.  Price. 
"Acoustic  Attenuation  Experiments  on  Subscale,  Cold-Flow  Rocket  Motors," 
AMER  INST  AERONAUT  ASTRONAUT  J,  Vol.  5,  No.  2 (February  1967),  pp.  272- 
80. 
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grid  which  was  used  in  the  problem  and  isometric  views  of  the  grid  (one 
for  each  acoustic  mode)  which  show  the  acoustic  pressure  distribution  in 
a vectorial  manner  for  each  standing  wave  solution. 

The  tabular  pressure  distributions  are  mainly  of  use  only  when 
detailed  quantitative  information  is  required  of  the  acoustic  pressure 
and  that  information  is  best  relegated  to  an  appendix.  The  isometric 
graphical  output  showing  the  acoustic  pressure  distribtuion  is  quite 
useful  for  a quick,  qualitative  view  of  the  nature  of  the  acoustic  wave 
structure  and  extensive  use  of  isometric  graphics  is  made  in  describing 
the  results  in  this  report. 

The  3D  NASTRAN  acoustics  program  results  are  structured  around  the 
order  of  tangential  solutions.  The  set  of  solutions  for  which  the  tan- 
gential order  is  zero  contains  all  pure  axial,  pure  radial,  and  combina- 
tion axial-radial  waves.  The  highest  mode  number  allowed  in  the  analy- 
sis was  normally  set  for  a value  of  ten.  A mode  number  in  excess  of  20 
would  be  needed  for  the  lowest  possible  radial  wave  solution  to  be 
reached.  Therefore,  all  zero  order  tangential  solutions  which  were 
obtained  were  of  axial  waves  only. 

Tangential  solutions  of  order  unity  include  pure  first  tangential 
waves,  combination  first  tangential-axial  waves,  and  combination  first 
tangential-axial-radial  waves.  Only  pure  first  tangential  and  combina- 
tion first  tangential-axial  wave  solutions  were  obtained  as  the  number 
of  modes  allowed  was  not  high  enough  to  permit  radial  solutions  to  be 
obtained.  Similarly,  second  order  tangential  solutions  include  pure 
second  tangential  waves  and  combinations  of  second  tangential,  axial, 
and  radial  waves.  As  with  the  sets  of  zero  and  first  order  tangential 
solutions  a mode  limit  of  ten  was  imposed  and  no  solutions  containing 
radial  wave  motions  were  obtained.  No  third  or  higher  order  tangential 
solutions  were  run. 

Isometric  graphic  displays  of  the  acoustic  pressure  distributions 
for  the  four  lowest  axial  frequencies  are  shown  in  Figure  4.  These  were 
obtained  with  the  assumption  that  the  nozzle  throat  is  closed.  The 
vertical  lines  in  the  figure  represent  the  relative  magnitude  of  the 
acoustic  pressure  at  each  grid  point.  The  pressure  distributions  shown 
are  the  acoustic  perturbations  about  the  mean  chamber  pressure.  The 
perturbed  values  are  presented  as  if  frozen  at  a point  in  time  when  the 
magnitude  of  the  maximum  acoustic  pressure  in  the  cavity  has  reached  an 
arbitrary  value  of  unity.  The  distribtuions  shown  in  Figure  4,  and  in 
all  similar  figures  in  this  report,  are  for  0 cm  web  burn. 

Acoustic  pressure  distributions  for  the  four  lowest  axial  modes 
with  an  open  throat  appear  in  Figure  5.  Two  notable  differences  between 
closed  and  open  throat  solutions  are  that  for  the  same  mode  number  the 
closed  throat  frequencies  are  higher  and  there  are  pressure  antinodes 
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15.26  Hz 


31.09  Hz 


47.80  Hz 


64.74  Hz 


FIGURE  4.  Acoustic  Pressure  Distribution  for  the 
Four  Lowest  Axial  Frequencies  - Closed  Throat. 

(0  web  burn.) 
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present  in  the  nozzle  end  of  the  motor  when  a closed  throat  is  assumed 
but  not  when  the  throat  is  assumed  to  be  open.  A complete  listing  of 
axial  frequencies  for  the  three  web  burn  distances  and  with  closed  and 
open  throat  assumptions  appears  in  Table  1. 


TABLE  1.  SRM  Computer-Predicted  Frequencies  (NASTRAN) 
Axial  Solutions  - Hz. 


I 


Axial 

Nozzle 

Distance  of  web  burned 

mode 

No. 

l nr  os u 

condition 

0 cm  (0  in. ) 

48  cm- (19  in. ) 

86  cm  (34  in.) 

1 

Closed 

15.26 

13.98 

16.19 

X 

Open 

12.06 

11.49 

11.32 

Closed 

31.09 

29.89 

32.19 

Open 

22.48 

23.92 

24.11 

■j 

Closed 

47.80 

47.03 

48.64 

J 

Open 

40.29 

38.66 

38.61 

A 

Closed 

64.74 

64.94 

65.72 

H 

Open 

57.62 

55.85 

54.26 

Experience  with  axial  mode  instability  in  solid  propellant  rocket 
motors  indicates  that  the  strongest  mode  is  normally  the  fundamental 
(first)  mode.  Since  the  acoustic  pressure  distribution  in  the  SRM  of 
the  first  axial  mode  is  of  interest  to  structural  engineers  and  to  those 
Interested  in  minimizing  POGO  effect,  first  axial  mode  acoustic  pres- 
sures are  tabulated  for  each  of  the  three  web  burns  and  for  both  nozzle 
throat  conditions:  closed  and  open.  These  data  are  presented  in 
Appendix  E. 

Graphic  acoustic  pressure  distribution  for  the  four  lowest  fre- 
quencies obtained  for  solutions  of  tangential  order  unity  are  shown 
using  a closed  throat  assumption  in  Figure  6 and  for  an  open  throat 
assumption  in  Figure  7.  It  is  characteristic  of  this  class  of  solutions 
that  acoustic  wave  activity  in  the  lower  mode  numbers  occurs  either 
primarily  in  the  slotted  portion  of  the  motor  at  the  forward  end  or  in 
the  annular  space  that  surrounds  the  nozzle.  It  is  also  characteristic 
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FIGURE  7.  Acoustic  Pressure  Distributions  for  the  Four  Lowest 
First  Tangential  Frequencies  - Open  Throat  (0  Web  Burn). 
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that  at  higher  mode  numbers  than  shown  in  this  report  the  acoustic  wave 
structure  also  involves  the  circular  bore  portion  of  the  motor  cavity. 

A complete  listing  of  frequencies  for  tangential  solutions  of  order 
unity  are  shown  for  the  three  web  burn  distances  and  for  closed  and  open 
throat  in  Table  2. 


TABLE  2.  SRM  Computer-Predicted  Transverse  Frequencies  (NASTRAN) 
First  Tangential  Solutions  - Hz 


Mode 

number 

Nozzle 

throat 

condition 

Distance  of  web  burned 

0 cm  (0  in.) 

48  cm  (19  in.) 

86  cm  (34  in.) 

1 

Closed 

201.6 

186.3 

178.0 

Open 

201.6 

NDa 

ND 

2 

Closed 

246.8 

189.1 

186.0 

Open 

253.6 

ND 

ND 

3 

Closed 

255.0 

224.1 

195.3 

Open 

255.0 

ND 

ND 

4 

Closed 

331.4 

255.9 

198.9 

Open 

334.6 

ND 

ND 

a 

Not  determined. 


Tangential  solutions  of  order  two  are  shown  for  the  four  lowest 
frequencies  with  a closed  throat  in  Figure  8 and  with  an  open  throat  in 
Figure  9.  As  in  the  case  of  the  first  order  tangential  solutions, 
acoustic  activity  of  second  tangential  waves  in  the  lower  mode  numbers 
is  primarily  confined  to  the  slotted  head  end  region  or  to  the  annular 
apace  around  the  nozzle.  Again,  as  with  waves  of  the  first  tangential 
class,  second  order  tangential  waves  also  occur  in  the  central  circular 
bore  of  the  motor  but  at  frequencies  higher  than  are  shown  here.  Fre- 
quencies for  tangential  solutions  of  order  two  are  shwon  for  the  three 
, burn  distances  in  Table  3. 
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FIGURE  9.  Acoustic  Pressure  Distributions  for  the 
Four  Lowest  Second  Tangential  Frequencies  - Open 
Throat  (0  Web  Burn). 
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TABLE  3.  SRM  Computer-Predicted  Transverse  Frequencies  (NASTRAN) 
Second  Tangential  Solutions  - Hz 


Mode 

Nozzle 

Distance  of  web  burned 

number 

tnroat 

condition 

0 cm  (0  in. ) 

48  cm  (19  in.) 

86  cm  (34  in.) 

1 

Closed 

218.4 

299.7 

291.2 

Open 

218.4 

NDa 

ND 

2 

Closed 

273.3 

304.4 

304.9 

Open 

273.2 

ND 

ND 

3 

Closed 

356.6 

333.1 

319.6 

Open 

356.6 

ND 

ND 

A 

Closed 

357.6 

391.8 

325.4 

Open 

358.8 

ND 

ND 

flNot  determined. 


First  and  second  tangential  results  for  48-  and  86-cm  web  burns 
were  obtained  only  for  the  closed  throat  condition  since  the  influence 
of  an  open  throat  diminishes  with  web  burn  and  the  cost  of  additional 
computer  time  did  not  appear  to  be  warranted. 


DISCUSSION  AND  CONCLUSIONS 


Initial  hand  calculations  of  SRM  acoustic  characteristics  have  been 
supplanted  by  more  accurate  3D  NASTRAN  results.  The  hand  calculations 
were  limited  by  assuming  that  the  nozzle  throat  was  an  acoustic  reflect- 
ing surface  (closed)  and  by  approximating  the  actual  motor  geometry  with 
an  equivalent  right  circular  cylinder;  the  main  advantages  of  hand 
calculations  were  speed  and  economy.  The  3D  NASTRAN  method  for  obtain- 
ing acoustic  characteristics  accounts  for  the  acoustic  effects  of  slots, 
tapers,  and  other  geometric  complications  within  the  motor  and,  thereby, 
generally  provides  more  accurate  acoustic  solutions  than  can  be  obtained 
by  simple  hand  calculations.  The  NASTRAN  technique  allows  the  option  of 
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treating  the  nozzle  throat  as  a closed  or  an  open  area  in  regard  to 
acoustic  wave  reflections.  Finally,  the  NASTRAN  program  also  provides 
detailed  data  regarding  acoustic  pressure  distribution  necessary  for 
performing  combustion  stability  analysis  of  the  motor. 

In  comparing  axial  acoustic  wave  frequencies  and  wave  structure, 
the  differences  between  the  closed  and  open  nozzle  throat  conditions  are 
apparent:  when  the  closed  throat  assumption  is  applied,  acoustic  waves 
impinging  on  the  nozzle  throat  are  reflected  and  a pressure  antinode  can 
exist  in  the  aft  end  of  the  motor;  use  of  the  open  throat  assumption 
results  in  a condition  in  which  a pressure  antinode  cannot  exist  in  the 
immediate  vicinity  of  the  throat.  The  effect  of  the  throat  is  greatest 
early  in  burn  when  J is  largest.  The  effect  of  the  open  throat  on  axial 
mode  frequencies  is  to  cause  a lower  frequency  to  exist  for  a given  mode 
number  than  for  the  closed  throat  frequency.  However,  the  influence  of 
throat  condition  on  frequency  diminishes  as  J decreases  and  as  mode 
number  increases. 

It  is  not  possible  at  this  point  to  know  which  set  of  axial  wave 
solutions  will  prove  the  most  accurate  for  the  SRM.  Therefore,  solu- 
tions for  both  throat  conditions  are  provided. 

Detailed  information  regarding  the  finite  element  grids  and  the 
fundamental  axial  wave  pressure  distribution  are  included  in  Appendices 
B,  C,  and  E for  use  by  the  reader  in  programs  which  are  involved  with 
the  effect  of  an  internal  acoustic  wave  on  the  motor  and  Shuttle  struc- 
tures. 

% 

The  tangential  wave  solutions  are  of  somewhat  more  complicated 
structure  than  the  axial  waves.  Early  in  burn  when  the  slotted  portion 
in  the  forward  end  of  the  motor,  the  cylindrical  centerbore,  and  the 
annular  space  surrounding  the  nozzle  are  geometrically  most  distinct 
from  each  other,  the  transverse  solutions  show  acoustic  activity  primar- 
ily in  the  extreme  forward  and  aft  portions  of  the  motor  for  the  lower 
frequencies.  As  the  frequency  increases,  transverse  waves  also  involve 
the  circular  bore  portion  of  the  motor.  The  tangential  solutions  show 
that  at  later  stages  of  burn  as  the  head,  center,  and  aft  portions  of 
the  motor  become  geometrically  less  distinct  from  each  other,  the  tan- 
gential wave  solutions  tend  to  couple  more  readily  between  one  portion 
of  the  motor  and  another. 

In  regard  to  the  POGO  effect,  with  its  approximately  50  Hz  upper 
frequency  limit,  the  only  acoustic  waves  in  the  SRM  likely  to  interact 
with  POGO  are  the  lowest  axial  waves.  No  transverse  wave  in  the  SRM  has 
a predicted  frequency  below  approximately  180  Hz,  therefore  waves  of 
that  class  are  outside  the  range  of  interest  to  POGO. 
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Acoustic  wave  pressure  distributions  for  all  inodes  run  on  the  3D 
NASTRAN  program  will  be  kept  on  file  at  NWC.  Should  a need  for  such 
Information  arise,  it  can  be  obtained  on  request  by  contacting  either  of 
the  authors. 

No  2D  analyses  of  the  SRM  have  been  made  to  date  as  a requirement 
for  2D  data  has  not  been  established.  However,  this  program  is  opera- 
tional and  could  be  used  to  obtain  higher  resolution  of  transverse 
acoustic  wave  characteristics  than  have  been  obtained  with  the  3D  pro- 
gram. 
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Appendix  A 

USE  OF  CLASSICAL  ACOUSTICS  OF  A RICHT  CIRCULAR  CYLINDER 
TO  ESTIMATE  ROCKET  MOTOR  FREQUENCIES 


Resorting  to  simple  methods  for  predicting  rocket  motor  acoustics 
which  involve  use  of  classical  acoustics,  simplifying  assumptions  regard- 
ing the  interior  geometry  of  the  motor,  and  application  of  a simple 
closed-form  algebraic  relation  which  allows  motor  frequencies  to  be 
calculated  quickly  by  hand  might  seem  antiquated  and  out  of  place  when 
compared  with  the  elegant  finite-element  methods  currently  available 
which  permit  the  frequency  and  acoustic  pressure  distribution  to  be 
calculated  to  virtually  any  desired  degree  of  precision.  However,  the 
simple  classical  approach  has  its  place  when  time  and  cost  are  at  a 
premium  and  when  approximate  estimates  of  acoustic  wave  frequencies  are 
sufficiently  accurate  at  least  on  an  interim  basis. 


The  first  acoustic  frequency  calculations  to  be  executed  at  NWC 
concerning  the  SRM  were  based  on  the  classical  acoustics  model  described 
below  and  the  assumption  that  the  actual  motor  geometry  can  be  described 
in  terms  of  an  equivalent  right  circular  cylinder.  The  results  were 
distributed  to  participants  at  early  meetings  concerned  with  assessment 
of  SRM  combustion  stability.  Since  the  hand  calculations  were  used  in 
early  discussions  of  SRM  combustion  stability,  they  will  be  discussed  in 
more  detail  than  has  been  done  previously. 


Acoustic  oscillations  in  a fluid  medium  are  pressure  oscillations 
of  small  amplitude  and  are  described  mathematically  by  the  classical 
wave  equation.  For  a right  cylindrical  cavity  with  closed  ends  and 
ideally  rigid  walls  the  acoustic  pressure  variation  can  be  calculated 
using: 


»,n,n 


tn,  n , n 


[AjCos(m$  - wt  - 6^)  + A^cosCmifi  - wt  - 


(A-l) 


$2)]  (Eq.  3 of 
footnote  7) 


in  which 


is  the  difference  between  local  and  space  averaged  pressure 
at  any  point  in  space  and  time 

Are  the  cylindrical  coordinates  with  the  origin  at  the 
center  of  one  end  of  the  cavity 


-7R.  D.  Smith  and  D.  F.  Sprer>ger.  "Combustion  Instability  on  Solid 
Propellant  Rockets",  Fourth  Symposium  on  Combustion,  Williams  & Wilkins 
Co.,  Baltimore,  19S3. 
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m,n,n 


A1,A2 

V62 


Radius  and  length  of  the  cavity 

Wave  numbers  characterizing  any  particular  mode  of  oscillation 

Bessel  function  of  order  m 

nth  root  of  the  equation  ~~  J (*x)  * 0 

dx  m 

(Some  values  are  given  in  Table  A-l) 

Arbitrary  independent  amplitude  constants 
Arbitrary  independent  phase  constants 


u Circular  frequency 

Every  possible  acoustic  mode  has  its  frequency  which,  for  a cylin- 
drical cavity,  can  be  calculated  using  the  following  equation: 


where  c is  the  velocity  of  sound  of  the  gas  in  the  cavity. 

% 

Any  particular  mode  of  oscillation  is  identified  by  the  wave  number 
in  each  of  the  three  directions,  axial  (nz) , radial  (n) , and  tangential 
(m).  Values  of  for  wave  numbers  up  to  3 are  given  in  Table  A-l. 
Where  only  one  wave  number  is  not  zero,  the  corresponding  mode  is  a pure 
mode.  For  example,  axial  acoustic  waves  have  n2  ^ 0,  m = 0,  and  n = 0. 
The  axial  wave  number,  nz,  is  expressed  as  a positive  integer.  Thus 
axial  mode  frequencies  are  given  by 


where  n ■ 1,  2,  3,  . . . 
z 

Likewise,  a pure  tangential  wave  frequency  is  given  by  the  relation 


where  is  determined  for  n ■ 0 and  m any  positive  integer.  In  the  case 

of  the  first  tangential  wave,  for  example,  a - 0.586. 

. mn 
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Tangential 
wave  No. . 


TABLE  A-l.  Values  of  a 


Radial  wave  number,  n 


2 


0.000  1.220  2.233 
0.586  1.697  2.714 
0.972  2.135  3.173 
1.337  2.551  3.611 


Although  all  combinations  of  pure  and  mixed  waves  are  possible,  it 
has  been  the  practice  to  deal  primarily  with  the  lowest  three  or  four 
frequencies  of  the  pure  modes.  The  results  for  the  SRM,  using  a speed 
of  sound  of  9.9  x 10^  cm/s  (3,250  ft/s),  are  shown  in  Table  A-2. 
Dimensions  used  in  the  calculations  are  shown  in  Table  A-3. 


TABLE  A-2.  NASA  SRM  Acoustic  Calculations 
(Preliminary  Hand  Calculations) 


Mode  No 


Acoustic  frequencies,  Hz 


AVJ  - 1 

Tangential 

Radial 

AXldX 

At  ignition 

At  burnout 

At  ignition 

At  burnout 

15.5 

363 

160 

755 

333 

31.0 

602 

265 

1,382 

609 

46.5 

828 

365 

• a • 

• • • 

62.0 

• • • 

• • • 

• • • 

• • • 
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TABLE  A- 3 


Dimensions  used  In  calculations: 


Motor  length  (interior) 


At  Ignition 


3,193  cm 
(104.75  ft) 


At  burnout 


3,193  cm 
(104.75  ft) 


Circular  perforation 
diameter 


160  cm 
(63  in.) 


363  cm 
(143  in.) 


NWC  TM  2822 


Append  I x 8 
FINITE  ELEMENT  GRIDS 


- 

> *• 

?*!FS 

- ■ •*—  • 

fcircrxirrirrr 

\»  r» 

zz: 

* Ui*  C 

• 

)•  M»  - 

rWt 

i«t 

t •« 

IN 

Forward  Portion 


I II IITI L L LU  UJJJJJXCir 

Center  Portion 


TTTTTTTT 

•»  — 

•m  — 

> 

ZJi 

llif  l||>l  MM  III*  III*  III* 

u- 

\s\ 

VlLruiuBtv.. 


Aft  Portion 


Finite  Element  Grid  for  0-cm  Web  Burn  Acoustic  Analysis 
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Aft  Portion 


Finite  Element  Grid  for  48-cm  Web  Burn  Acoustic  Analysis 


Aft  Portion 


Finite  Element  Grid  for  86-cm  Web  Burn  Acoustic  Analysis 
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Appendix  C 

FINITE  ELEMENT  GRID  COORDINATE  DATA 


The  following  tabulations  provide  all  the  necessary  data  for  deter- 
mining the  finite  element  grids  used  in  the  SRM  acoustic  calculations. 
Three  tables  are  presented,  one  for  each  web  burn. 

The  left  hand  column  identifies  the  order  of  the  card  in  the 
sequence. 

The  AXSLOT  card  contains  the  gas  density  and  bulk  modulus,  the 
tangential  number  and  two  fields  that  have  a default  slot  width  and 
number  of  slots.  The  SLBDY  card(s)  lists  the  grid  points  along  the 
slot-gas  cavity  border. 

The  majority  of  cards  are  of  two  types:  element  cards  and  grid 
cards.  These  are  described  in  the  following: 

CAXIF2  are  centerbore  elements  along  the  centerline. 

CAXIF3  are  three-sided  fluid  elements 

CAXIF4  are  four-sided  fluid  elements 

CSL0T3  are  three-sided  fluid  elements  in  the  slots 

CSL0T4  are  four-sided  fluid  elements  in  the  slots 

GRIDF  points  form  the  corners  of  the  CAXIF  elements 

GRIDS  points  form  the  corner  of  the  CSLOT  elements 

The  second  column  of  the  element  cards  contains  the  element  identi- 
fication number.  The  next  two,  three,  or  four  columns  contain  the  grid 
identification  numbers  of  the  corner  points  of  that  element. 

The  second  column  of  the  grid  cards  is  the  grid  identification 
number.  The  third  column  is  the  distance  from  the  centerline  (R)  in 
inches.  The  fourth  column  is  the  axial  distance  from  the  reference 
point  (z)  in  inches. 
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TABLE  C-l 
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Grid  Coordinate  Data  for  0-cm  Web  Burn 
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Appendix  D 

EFFECT  OF  J ON  ROCKET  MOTOR  FREQUENCY 


A series  of  tests  involving  cold  gas  flow  through  a small  scale 
rocket  model  was  conducted  to  determine  the  effect  of  the  nozzle  on 
acoustic  losses.*  The  model  was  excited  at  its  fundamental  axial  fre- 
quency by  an  acoustic  driver.  Although  there  was  considerable  scatter 
in  some  of  the  data,  it  was  noted  that  the  frequency  generally  decreased 
as  the  ratio  of  the  nozzle  throat  area  to  motor  gas  channel  area  (J) 
increased.  Figure  D-l  is  based  on  data  in  footnote.*  In  this  figure 
the  data  has  been  fitted  with  a straight  line  using  a least  squares 
technique.  The  frequency  scale  has  been  normalized,  using  a frequency 
of  5A5  Hz  as  the  normalizing  factor.  The  boundaries  of  the  figure  have 
been  expanded  to  show  a range  of  J from  0 to  1 and  the  fitted  straight 
line  has  been  extrapolated  to  higher  values  of  J than  were  used  in 
obtaining  the  experimental  data. 

While  the  experimental  data  seems  to  extrapolate  to  the  approximate 
vicinity  of  the  classical  frequency  prediction  for  J*l,  it  should  be 
noted  that  the  effects  of  flow  are  entirely  absent  in  the  classical 
prediction  and  that  the  experimental  data  is  for  values  of  J for  which 
the  gas  velocity  (Mach  number)  is  relatively  small.  Thus,  the  existing 
experimental  data  does  not  appear  to  provide  information  on  the  effect 
of  gas  flow  on  acoustics. 

The  results  suggest  that  the  fundamental  frequency  of  a rocket 
motor  is  a function  of  J.  Therefore,  the  assumption  that  the  nozzle 
throat  has  no  direct  influence  on  an  axial  acoustic  wave  may  be  in 
error. 


*The  test  results  conducted  by  Buffum,  et.al,  were  reported  in  1967 
in  the  AIAA  Journal  (see  Ref.  6 on  page  9 of  this  memorandum). 


NORMALIZED  FREQUENCY  - f/f.  (f.  - 545  Hz) 
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Appendix  E 

RELATIVE  ACOUSTIC  PRESSURE  AT  GRID  POINTS 
FOR  FIRST  AXIAL  MODE 

The  relative  pressure  at  grid  points  1 through  6 are  listed  across 
row  1.  Row  2 contains  the  pressures  for  grid  points  7 through  12,  etc. 

If  a number  was  not  used  in  the  finite  element  grid,  the  corresponding  point 
in  the  pressure  listing  .will  be  blank. 

T Headings,  Column  Designations 


T1 

T2 

T3 

R1 

R2 

R3 


N 

N+l 

N+2 

N+3 

N+4 

N+5 


N » number  in  left  column 

The  type  (S)  column  is  not  applic~ 
able  in  this  tabulation. 
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